Weyl problem and Casimir effects in spherical shell geometry 
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We compute the generic mode sum that quantifies the effect on the spectrum of a harmonic field 
when a spherical shell is inserted into vacuum. This encompasses a variety of problems including 
the Weyl spectral problem and the Casimir effect of quantum electrodynamics. This allows us to 
resolve several long-standing controversies regarding the question of universality of the Casimir self- 
energy; the resolution comes naturally through the connection to the Weyl problem. Specifically we 
demonstrate that in the case of a scalar field obeying Dirichlet or Neumann boundary conditions on 
the shell surface the Casimir self-energy is cutoff-dependent while in the case of the electromagnetic 
field perturbed by a conductive shell the Casimir self-energy is universal. We additionally show that 
an analog non-relativistic Casimir effect due to zero-point magnons takes place when a non-magnetic 
spherical shell is inserted inside a bulk ferromagnet. 

PACS numbers: 03.70.+k, ll.10.-z, ll.10.Gh, 42.50.Pq 
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I. INTRODUCTION 



Many physics problems require the evaluation of the 
sum 

/>oo 

Z{Q) = Y j F(q a ) a / F{q)G(q)dq (1) 
•'0 

where q a is the wavevector spectrum for a harmonic field 
confined to a region and F(q) is a function that for q 3> Q 
goes to zero monotonically and fast enough that the mode 
sum converges (the essence of the second step in (JTJ will 
be explained shortly). The parameter Q is the cutoff 
wavevector whose physical meaning depends on the sys- 
tem in question. When the function F is the contribution 
of a mode to a thermodynamic property of a harmonic 
field of wavevector q such as the entropy, energy, or free 
energy, the corresponding sums (JTJ describe the thermo- 
dynamics of black-body radiation in a cavity or,equiva- 
lently, the thermodynamics of a harmonic solid [l| . Here 
the Planck distribution enters the function F(q) guaran- 
teeing the convergence of the sum ([1]) and the tempera- 
ture T serves as the wavevector cutoff Q. For a macro- 
scopic region of typical size a the condition T > 1/a is 
satisfied and there are many terms in the 1/a < q a < T 
range contributing to the mode sum. 

The function F can also represent the zero-point en- 
ergy of a harmonic field with linear dispersion law; how- 
ever it is useful to amend the dependence F = q/2 at 
large wavevectors so that F goes to zero at q large. This 
assures that the sum (JTJ) converges, and helps classify and 
resolve the divergences at high wavevectors that might 
otherwise occur; physically this means that a material 
boundary surrounding the region in question becomes in- 
visible to short wavelength radiation [3 Q and thus does 
not perturb the spectrum at highest wavevectors. In this 
example the cutoff wavevector Q is a property of the ma- 
terial of which the boundary is constructed. Mode sums 
like (JTJ are also encountered in Fermi systems [lj. 

In a macroscopic system the discrete spectrum can be 
approximated by a continuous one and characterized by 



a density of states (DOS) G(q) so that the mode sum can 
be represented as an integral [l| as shown in the second 
step in (JTJ. The function G(q), to be referred to as the 
Weyl DOS, represents the continuum approximation to 
the exact DOS Q{q) = J2 q ^(° ~~ ?<*)) ^ plays a central 
role in the physics of finite-sized systems. The significant 
feature of this replacement is that G is independent of 
F: it encodes the physics of the field and the boundary 
condition, while F is the specific aspect of the system 
that is being studied. 

The large-g behavior of the Weyl DOS G(q) can be 
inferred from the behavior of Z(Q) for large cutoff by in- 
verting the integral transform ([TJ. These leading terms, 
known as the Weyl series [J], have geometric interpre- 
tation. The Weyl series is the origin of the divergences 
that occur in attempts to calculate the Casimir effect 
[3 Q without use of a cutoff; such divergences can repre- 
sent real effects of the presence of the physical cutoff Q 
ii- 

Since the Weyl DOS G(q) is independent of the cutoff 
procedure, it can be extracted from the mode sum calcu- 
lated using any convenient function F(q). Kac [l(j and 
Stewartson and Waechter [ll[ demonstrated the utility 
of the Gaussian function F = cxp(—q 2 t) (Q = 1/vt) 
and applied it to the case of a two-dimensional region. 
However the exponential F = exp(-qt) and power-law 
F oc q~ s choices are just as good. The latter is employed 
in the zeta function regularization method [13]; there is 
no cutoff scale Q, and one studies instead the role played 
by the parameter s. 

The goal of this paper is to demonstrate how to com- 
pute the mode sum for generic F(q) and thus to infer the 
change in the Weyl DOS due to the insertion of a three- 
dimensional spherical shell into vacuum. Along the way 
we will also be able to resolve long-standing controver- 
sies regarding the role of short- wavelength cutoff on the 
Casimir self-energy. 
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II. ONE-DIMENSIONAL GEOMETRY 

To illustrate the concepts we begin with an example 
of a one-dimensional interval of length a and a scalar 
harmonic field satisfying the Dirichlet boundary condi- 
tions at the interval ends. The spectrum is given by 
q n = irn/a,n — 1,2,.... Then employing the Euler- 
Maclaurin summation formula flij the mode sum (JlJ can 
be transformed into 



Z(Q,a) = Y,F(—) = - / F(l)dq- 
a nJ ° 



F(0) 



ttF'(O) tt 3 F( 3 )(0) 



12a 



720a 3 



■(2) 



assuming that F(y) is not singular at the origin and that 
it and all of its derivatives vanish at infinity. In view of 
Eq. ([l} this is consistent with the Weyl expansion of the 
form 



n , \ a %) 

7T Z 



w6'(q) w 3 S^(q) 



12a 



720a 3 



(3) 



with the understanding that the delta-function and its 
derivatives are only used as a shorthand to indicate 
that the mode sum is sensitive to aspects of the long- 
wavelength part of F(q). 

In this example the first two terms of the Weyl ex- 
pansion §3§ have geometrical relevance: the leading term 
is proportional to the one-dimensional "volume" a; the 
corresponding term of the mode sum ([2]) is the only one 
that requires a cutoff. The next order delta-function term 
represents the effect of the edges of the interval. Our ex- 
pression Q contains more terms than usually kept [J; 
typically one keeps only the geometric terms and deals 
with the rest separately. Our attitude here is that the 
DOS can have as many terms as needed as long as for 
any physical quantity the outcome can be presented in 
the integral form like in Eq.([T]). 

The coefficient of the second term in ([2]) (which gives 
rise to the delta- function term in ([3])) has a special place 
within the theory because it does not depend on the 
wavevector cutoff Q or on the macroscopic length scale 
a; it is just due to the presence of the boundary. We will 
call this coefficient the Kac number /C; apparently this 
number was first computed (in any context) by Kac (lpj 

who found that IC d J 2 = 1/6 for a simply-connected two- 
dimensional region enclosed by a smooth Dirichlet curve. 
For the one-dimensional Dirichlet interval Eqs.® and §S§ 

imply that IC d=1 = —1/2. The Kac number gives geomet- 
rical information about the boundary and its topology. 
The Kac term does not contribute into the zero-point en- 
ergy, but it has other measurable consequences because 
K, reflects the change in the number of states due to intro- 
duction of a boundary. Then the classical equipartition 
theorem [l[ implies that the energy of a region contains a 
universal KT term [TEl ]; in one dimension the —T/2 piece 
in the energy is the leading finite-size contribution. 



The derivative terms in ([2]) and (j3]) also have a special 
place within the theory because they only depend on the 
macroscopic length scale a and do not depend on the 
cutoff. They are responsible for the Casimir effect and 
its generalizations. The generalized Casimir effect will 
be defined as a change of the value of the mode sum as 
a result of the introduction of a boundary or boundaries. 
For example, the change of the mode sum as a result of 
inserting into vacuum of two Dirichlet points separated 
by a distance a will be given by 



AZ(Q,a) = -F(0) 



ttF'(O) tt 3 F( 3 )(0) 



12a 
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(4) 



because the mode sum for the vacuum contains only the 
leading term of ([2]); the first term is due to the insertion 
of two Dirichlet points. Choosing F(q — > 0) — > q/2 we 
recover the well-known — 7r/24a Casimir attraction [l6[ 
given by the second term of We also note that if 
there were a definite function F(q) that we were study- 
ing which had a nonzero third derivative at zero wavevec- 
tor, Eq.Q would describe the consequences which would 
be an example of generalized Casimir effect. Since the 
mode sum ((4]) only contains odd derivatives of F, the 
generalized Casimir effect is absent for any function that 
for small wave vectors vanishes as an even power of the 
wavevector. Specifically, this rules out the possibility of 
the Casimir effect with " non-relativistic" dispersion law 
F(q — > 0) ex q 2 ; the same conclusion holds in the parallel 
plane geometry in three spatial dimensions [17j . How- 
ever below we will demonstrate a possibility of a non- 
relativistic Casimir effect in spherical shell geometry. 

The mode sum @ is consistent with the Weyl expan- 
sion of the form 



AG(q) = -S(q) 



n5'(q) 7r 3 (5( 3 )(g) 



12a 



720a 3 



(5) 



Since the original Dirichlet interval can be viewed as a 
one-dimensional sphere while the configuration with two 
Dirichlet points represents a one-dimensional spherical 
shell (both of radius a/2) it is instructive to compare the 
expression for the DOS (|5|) accumulating the effect of 
the field modes both inside and outside of the Dirichlet 
interval with that given by Eq.© which only includes 
the effect of the interior modes. We then observe that 
the geometric part of the DOS given by the first two 
terms of © and by the first term of (O is additive: the 
bulk a i 'tt term present in the original DOS © is can- 
celled between the interior and exterior modes while the 
effect of the edges given by the delta-function term in 
(O is doubled. Second, the derivative terms of the DOS 
expressions © and (|5|) are identical which is an indica- 
tion of non-additivity of the generalized Casimir effect. 
These properties illuminated in the one-dimensional ex- 
ample hold more generally for the standard Casimir effect 
01- 
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III. THREE-DIMENSIONAL GEOMETRY 

We now proceed to a calculation of the counterparts 
of Eqs.(j4]) and (0 for a spherical shell in three dimen- 
sions starting with the scalar problem which is technically 
more involved than its electromagnetic counterpart. 



A. Scalar problem 

Special cases of the general mode sum have been con- 
sidered previously. 

* The sum fl} over the modes of a scalar field inside 
a Dirichlet sphere was first calculated by Waechter [i~8| . 
using the choice F(q) = exp(— q 2 i). The geometric parts 
of the DOS for a Dirichlet shell then can be extracted 
by use of the additive property of the geometric parts of 
the DOS. Below we will calculate the generic mode sum 
(HJ) and the DOS directly in the shell geometry and find 
some difference from Waechter's results. 

* The mode sum for the case F = q/2 corresponds to 
the zero-point energy of a scalar field. The change in the 
zero-point energy resulting from the imposition of a new 
boundary is the scalar Casimir effect. This was first stud- 
ied for the case of a spherical shell in three dimensions by 
Bender and Milton [ljj . The leading terms of the Weyl 
DOS give rise to ultraviolet divergences which were elimi- 
nated by means of a regularization technique. In order to 
study these divergences we will introduce a cutoff, so that 
F = (q/2)C(q/Q) (where C(q/Q) is small for large argu- 
ment). The divergences become cutoff-dependent terms; 
we will argue that these terms have physical meaning. 

To the best of our knowledge the connection between 
the Weyl problem and the Casimir problem in the case 
of a spherical shell has not been clearly established pre- 
viously. 



1. Contour integral representation of the mode sum 

In spherical geometry the modes are not spaced regu- 
larly in wavevector, and their determination entails solv- 
ing a transcendental equation. We can make use of this 
equation to represent the mode sum as a contour inte- 
gral. If q n and q p are the zeros and poles, respectively, of 
a function tp(g), the sum of the values of a function F(q) 
over these sets can be calculated by 



of a spherical shell of radius a that imposes Dirichlet 
boundary conditions. This will be written in the form 



1 f d 

— f F{q)—hx i p{q)dq = Y J F{qn)-Y. F ^v 



(6) 



where the contour C encloses the values being included 
in the sum; the function F(q) must be analytic inside 
the contour [2(1 HJ , and we assume that it is very small 
where the contour cuts the real q axis. 

We wish to calculate the change in the mode sum aris- 
ing from the introduction into a previously empty space 



AZ(Q,a) 



2=0 



2 v y(F(qi <k )-F(qj^)) 



k 



1=0 



vZ u (Q, a) 

_ (7) 
where qi t k and q^k are the spectra of the system with and 
without the shell, I labels the relevant spherical harmonic 
and the order of the corresponding Bessel function, and 
k labels the successive modes for given I. The factor 
2v = 21 + 1 accounts for mode degeneracy. 

To use the contour representation (J6j) to calculate the 
partial sum Z ll (Q, a) we construct the function ip(q) that 
vanishes when q is a mode wavevector of the modified 
system and has poles when q is a mode wavevector of the 
unperturbed system. To avoid dealing immediately with 
the complications of a continuous spectrum, consider the 
modes of a scalar field interior to a large sphere of radius 
b (with Dirichlet boundary conditions on the sphere), and 
the change in this spectrum caused by the introduction 
of a concentric Dirichlet shell of radius a < b. Then 

ip. v {q) = irqa[J v {qb)H£\qa) - J w ( 9 a)ff«(g&)]^M 

where J v (z) and (z) are the Bessel and Hankel func- 
tions, respectively The wavevectors of the modes 
interior to the inner sphere are the zeros of the factor 
J u (qa) 1 and similarly the modes of the unpartitioned 
sphere are the zeros of J u (qb) (and the poles of ip)\ 
the Bessel function combination vanishes when q is the 
wavevector for a mode confined to the region between 
the two spheres. The factor nqa has been introduced 
so that tp approaches unity for large imaginary q. This 
introduces a spurious mode at q = 0, which should be 
ignored. 

When q has an imaginary part and b >> a, J„(qb) is 
exponentially large and (qb) is exponentially small, 
so that tp reduces to ■naqJ v (qa)Hi, (qa). In this way 
the limit b — > oo can be taken. The net effect is the 
replacement [J v (qb)H^\qa) — J v (qa)H^ (qb)]/ J v (qb) — > 
H^iqa). 

Setting q = ivy /a, the partial sum Z v (Q,a) now has 
the form 



5sy>0 



F{ivy/a) d 
ri dy 



]n(2vyI v {vy)K v (vy))dy (9) 



The original integration contour C encloses the positive 
real q axis (but not including the origin, to exclude the 
spurious mode) in the counterclockwise direction. We 
can take it to enclose the whole right half-plane 5Rg > 0, 
which in the new variables is the upper half-plane 3y > 0. 
The part of the contour that lies at large ^sy makes no 
contribution because the cutoff function is very small 
there. This is similar to the starting point of others 
[22^ . but we have left this in the form of a contour in- 
tegral. We observe that the function F(q), which goes to 



4 



zero for large real q, becomes F{iyv / a) on the imaginary 
axis, and is not necessarily a "cutoff" function anymore; 
in particular, the exponential function has constant mag- 
nitude. For this reason it is important to keep the whole 
contour. 

We introduce the uniform asymptotic expansion of De- 
bye [H valid for v > 1: 



ln[2Kl + y 2 ) l ' 2 I v {vy)K v {vy)] = -^{-±- 

-i J 3 } + R{y, v) 
i + v ) > 



(i + y 2 ) 2 ( 



(10) 



The remainder function R(y, v) is of order v for large 
v and of order y~ A for large y. It gives a hnite con- 
tribution to ([7]) and thus may be disregarded while we 
consider the more problematic leading terms. With this 
approximation 



Fiivyja) d \ y 



<3y>0 

6 



dy L v/T + y 2 



— f — 

%v 2 1 1 + y 2 



: }]dy 



y 

(ii) 



(1 + 2, 2 ) 2 (1 + y 2 ) 3 . 

We observe that the only singularity of the integrand 
inside the contour is at y = i (that is, at q — vja in 
terms of the original variables); it is a superposition of 
poles of multiplicity 1,2, and 3. A little algebra leads to 



U a 

+ ^D-n-) + -^(-) + ^''(-)]+... 

64a z — ' a a a a a 

1=0 

(12) 

where the primes indicate differentiation with respect to 
the argument of F. 

Similar to the one-dimensional example studied above, 
the mode sum (Tj~2"|) can be understood with desired ac- 
curacy with the help of appropriate generalization of the 
Euler-Maclaurin summation formula [blj 



v 1 r°° i 



7 



1=0 



with the result 



5760 



/'"(0) 



(13) 



AZ(Q,a) = ~a 2 F(q)qdq- 
Jo 



F(0) 17F"(0) 



24 



1920a 2 



(14) 

where we restricted ourselves to the derivative terms 
whose order does not exceed the second. 



2. Analysis of the mode sum 

The result (|T4"]) is consistent with a change in the DOS 
of the form 



2 m , i7<y"(«) 



AG(q) = -a 2 q-^ + 



24 1920a 2 



(15) 



There is no bulk term proportional to a 3 in (fl4|) and (fl5|) 
because the volume of the system is not changed by the 
insertion of the shell. Even though approximations were 
made, the first two terms of (fT5j) are exact: the second 
of the sums in (IT2|) only contributes to the last term of 
(fT4| which also receives contribution from the first sum 
in (fl~2"|) . Consistent with Weyl's expansion, the leading 
terms of (|14p have geometrical interpretations. 

The first term is proportional to the area of the shell. 
Waechter 18], considering only the modes inside the 
shell, obtained a result half as large (in our evaluation, 
the shell perturbs the modes on both sides; the results 
are in agreement). This term has played a role in previ- 
ous treatments of the Casimir energy: when F(q) = q/2 
for all g, the mode sum (|12p becomes a divergent expres- 
sion. In previous publications [19l . |22j it was evaluated 
(to zero!) by employing the zeta regularization technique. 
We will discuss this approach below. 

The Kac number (the coefficient of S(q)) is K. d J 3 = 
— 1/24; it is negative because the modes of wavelength 
larger than a have been suppressed by the introduction of 
the shell. Modes inside and outside the shell are equally 
affected, so that a calculation that considers only the 
modes inside the sphere would give a Kac term that is 
half as large Waechter's calculation [lH overlooks 

the Kac number. 

We observe that the first derivative term is absent from 
Eqs.ffTJJ and (|15l) which is an indication that the cutoff- 
independent part of the Casimir effect has its origin in the 
remainder term R(y, v). Since this term does not require 
a cutoff, the existing treatment [l9[ is adequate and will 
not be repeated here. Thus the Casimir self-energy of 
the Dirichlet shell is given by 



t ~ 2 



q 2 C(q/Q)dq 



BM 



(16) 



where we wrote F(q) = (q/2)C(q/Q) with C(q/Q) rep- 
resenting the physical cutoff function determined by the 
transmission properties of the boundary and satisfying 
the conditions C(0) = 1 and C(oo) = 0. The l/a 
dependence of the cutoff-independent term of the self- 
energy (|16[) is dictated by dimensional analysis, and 
the numerical constant BA4 was computed by Bender 
and Milton [l9j]. The leading cutoff-dependent term 
of the self-energy has to be viewed as contributing 
(— l/87r) / °° q 2 C(q/Q)dq ~ — Q 3 into the bare coefficient 
of the surface tension of the shell, considered as material 
membrane. This would give rise to an outward stress on 
the sphere surface; without a cutoff, it would be an in- 
finite stress. This realization, originally due to Deutsch 
and Candelas 0, [H| , was recently re-expressed by Gra- 
ham and co-workers [25| and by Barton [26|. The im- 
plication is that in an experimental situation, a curved 
boundary might give rise to a cutoff-dependent contribu- 
tion to the physically measurable Casimir force (27| . 

Neumann boundary conditions for the scalar field can 
be discussed in the same way. The boundary condi- 
tion replaces I v (y) by ^fydj ' dy(^jyl v (y)) and K v {y) by 
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y/yd/ 'dy{-y/yK v (y)), and leads to an expression for the 
change in the density of states similar to Eq. (fT5|) 



respectively. For example, consider the evaluation of the 
sum (fT12j). With the choice F(q) = q~ s , this becomes 



AG(q) 



a 2 q 



7%) 976"(g) 
24 1920a 2 



(17) 



Since the effect of Neumann boundary conditions is to 
remove a constraint (continuity of the held across the 
shell) rather than to add one, every term has the opposite 
sign relative to the Dirichlet case. Similar to the Dirichlct 
case the Casimir self-energy is dominated by the cutoff- 
dependent contribution proportional to the area of the 
sphere. 



3. Non-relativistic Casimir effect 

Since Ecis- fT^]) and (|15l) contain second-order deriva- 
tive terms, they demonstrate the possibility of a non- 
relativistic Casimir effect. Since the non-relativistic 
Casimir effect of quantum electrodynamics is exponen- 
tially suppressed by the rest energy of elementary parti- 
cles, here we discuss its condensed matter analog. Even 
though the most interesting universal part of the ef- 
fect, as explained below, is small compared to its cutoff- 
dependent piece, the calculation is of interest as a proof 
of principle and because the overall effect is not small. 

Spin waves in ferromagnets are zero chemical poten- 
tial Bose excitations having a dispersion law whose long- 
wavelength limit is "non-relativistic" , u> = 70 2 |28[. Then 
the zero-point energy per mode is w/2 = -fq 2 /2. Let us 
now assume that there is a non-magnetic spherical sur- 
face of radius a (impenetrable to magnons) embedded 
inside bulk ferromagnet. The change in the zero-point 
energy due to the presence of the sphere can be found by 
substituting F(q) = (70 2 /2)C{q/Q) into Eq.(HU) 

(18) 

As in previously studied cases, the effect is dominated 
by the leading cutoff-dependent term reducing the bare 
coefficient of surface tension of the surface by an amount 
of the order 7Q 4 . The sub- leading term manifests itself 
as a universal finite-size correction with the 7/a 2 depen- 
dence dictated by dimensional analysis. The higher-order 
terms of the Debye expansion do not contribute anything 
further. 



4- Connection to the zeta function regularization method 

The zeta function regularization method evaluates the 
mode sum for the choice F(q) — q~ s . For sufficiently 
large s, the sums and integrals are convergent at large 
q; the result, referred to as the spectral zeta function, 
is analytically continued to physically relevant s. The 
relativistic and non-relativistic Casimir energies and the 
Kac number are the s = — l,s = — 2 and s — cases, 



AZ(s,a) 



-a s ( 



1 



a s (2 s - 1 -l)C( S -l) 
~' lS (5s + 6)(2 s+1 -l)C(s- 



G4 



1 



.(19) 



where £(2;, y) and (,{x) are the Hurwitz and Riemann 
zeta functions, respectively 
not convergent for s < 2, 



23] 



Although the sum is 
the zeta functions are de- 
fined by analytic continuation, and for physically inter- 
esting cases the spectral zeta function (TT91 reproduces 
all previously found universal results. Indeed, employ- 
ing C(-l) = -1/12, C(-2) = 0, and C(~3) = 1/120 
|23j we find AZ(-l.o) = (for the ordinary relativis- 
tic Casimir energy), AZ(-2,a) = 17/960a 2 (for the 7/2 
coefficient of the non-relativistic Casimir energy in Q18p). 
and AZ(0, a) — —1/24 (for the Kac number appearing 
in Eqs.CHD and ((15)1 ). 

The cutoff-dependent terms are also represented in the 
spectral zeta function. Eq. ((T9")) has a pole of residue 
— a 2 at s = 2 which corresponds to the leading cutoff- 
dependent term in the mode sum (|14l) . Indeed, evaluat- 
ing the latter with the choice F = q~ s and lower inte- 
gration limit q — e (because there is a lowest wavevector 
mode) we find — a 2 e 2 ~ s / (s — 2) which is a pole of residue 
—a 2 at s = 2 in agreement with the pole structure of (|19p . 
We claim that in general the cutoff-dependent terms in 
a mode sum are represented by the poles of the spectral 
zeta function that have to be bypassed on the way from 
s > 2 (or wherever the sums and integrals actually con- 
verge) to the physical case that is being studied (s = — 1 , 
s = —2, or s — in the above examples). Specifically, 
a pole of residue A at s — a signals the presence of an 
Aq a ~ x piece in the DOS which via ([T]) implies the pres- 
ence of cutoff-dependent contribution into the mode sum. 
We will insist that the cutoff-dependent terms have phys- 
ical meaning; they are a part of the Casimir energy (or 
whatever property is being studied) that is larger than 
the cutoff-independent contribution that is usually cal- 
culated. 

This connection explains why the zeta regularization 
fails for the calculation of the scalar Casimir effect due 
to a Dirichlet circle [l{|. The value of the Casimir term 
should be given by the case s = — 1 for the spectral zeta 
function AZ(s); however it has a pole there. The exis- 
tence of this pole is directly related to the logarithmically 
divergent AZ{Q) found by Sen [Ilj]. 

With this understanding, the zeta regularization is a 
powerful method for evaluating the mode sum. For ex- 
ample, suppose that (for the case of a spherical Dirichlct 
shell introduced into a scalar field) we wished to evaluate 
the change in the mode sum for the function F(q) = q 1 ^ 2 , 
which is the case s = —1/2 of the general problem 
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F(q) = q~ s . The singular behavior at q = is not con- 
sistent with the assumptions of the Euler-Maclaurin ex- 
pansion, so (fl~5|) does not apply. Since going from large s 
to the case of interest, s = —1/2, one encounters a pole 
at s = 2, there exists a corresponding cutoff dependent 
contribution. As above we will have to introduce a phys- 
ical cutoff function C. The value of (fT9 
cutoff-independent part. Thus for this example 



at s = — h is the 



AZ{Q,a) = -a 2 / q 3 / 2 C(q/Q)dq 



- 1/2 ((2- 3/2 -i)C(-5) + ^^ c 4 ) P ) 



B. Electromagnetic problem 

The most relevant mode sum problems involve find- 
ing the counterparts to Eqs. (fT4|) and (IT5t for the case of 
electromagnetic vacuum perturbed by a conductive shell. 



1 . Background 

The Weyl problem for the case of electromagnetic held 
in a cavity was considered by Balian and Bloch and by 
Balian and Duplantier [3(| ■ Boyer [3l[ demonstrated that 
the Casimir energy is cutoff-independent, varies as 1/a, 
and gives an outward stress on a sphere. Boyer's result 
has been re-derived in several complementary calcula- 
tions hi m, m . 

Unlike the case of the scalar held, the Casimir energy 
for the electromagnetic case is cutoff-independent. This 
is a special property of the vector character of the electro- 
magnetic field. Balian and Duplantier [l5[ argued that 
this is true for an arbitrarily shaped smooth conductive 
shell. Nearly all past work was criticized by Candelas [l3| 
who argued that the Casimir energy of a conductive shell 
is cutoff-dependent and that this cutoff-dependence can- 
not be explained in terms of the Weyl problem. Below, 
by computing the generic mode sum which encompasses 
both the Weyl and Casimir problems, we settle the con- 
troversy. 



Z (£M) = 



3y>0 



F{ivy/a) d 
iri dy 



(Ml-KW] 2 }) dy 



where 



°v(y) = yiv{y)K v (y) 



(22) 



(23) 



and the prime in (1221) indicates differentiation with re- 
spect to the argument of o v {y). If we choose F(q) = 
(q/2) exp(— g/Q) then Eas.(|2"T ]) -l|2"3" |) reduce to an expres- 
sion for the energy analyzed by Milton, DeRaadt and 
Schwinger [321 ] . 

The subsequent analysis mirrors the steps undertaken 
while treating the scalar version of the problem. In the 
present case the Debye expansion (TIT))) amounts to ap- 
proximation [32| 



KM? 



1 



4i/ 2 (l + y 2 f 

for v large. Evaluating the resulting contour integral 
Fiivyja) d ( I 



(24) 



Z (£M) _ _ 



3y>o 47r ^ 2 dy V (1 + 2/ 2 ) 3 



dy (25) 



and substituting the outcome into (f2"Tj) we find the elec- 
tromagnetic counterpart of Eq. (TT2"j) 



1 00 

AZ^(Q,a) = --L£[3F'(-) 



32a 

- 3- F"(~) + ^F'" (-)] + ... (26) 
a a a a 

We see that Eq. (|2"6"|) does not have an analog of the first 
sum in (|I2[) . indicating that a formal (Q = 00) treat- 
ment would face weaker divergences than in the scalar 
case. If we choose F(q) = q/2 for all q, then Eq.(|26|) 
reduces to the divergent expression for the Casimir en- 
ergy found by Nesterenko and Pirozhenko [22j where the 
divergence was evaluated by employing the zeta function 
regularization method; this time it gave a universal 1/a 
result. However, the evaluation is even simpler if F(q) 
contains an ultraviolet cutoff, because then the sum is 
always convergent. 

To use the Euler-Maclaurin summation formula (|13l) . 
write 



2. Analysis of the mode sum 

A derivation of an expression for the Casimir energy 
of the spherical shell beginning from the contour in- 
tegral representation (j6]) was given by Nesterenko and 
Pirozhenko (22j; its generalization to the case of generic 
mode sum requires only a few changes. Therefore we only 
quote counterparts of Eqs.® and ([9]) 



AZ( £ ^(Q,a)=J2uZ( £ ^(Q,a) 



(21) 



j>.]«5>.] -3F'(0) 



1=1 



1=0 



Then the mode sum ([26|) becomes 



F(0) 3F'(0) 
4 + 32a 



(27) 



(28) 



z=i 



The noteworthy features of this expression, compared to 
its scalar counterpart (fT4|) . are lack of the F"(0) term 
and the lack of dependence on the cutoff scale Q. The 
cutoff function did play a role, however: from the large v 



7 



and y dependences of (|24[) we can see that in the absence 
of the cutoff function, the integral (|2"21 converges and the 
sum (|21[) diverges. The effect of the cutoff function is to 
prevent the change in variables that would allow doing 
these calculations sequentially. 
The DOS is given by 

A G (-)(,)^ + ^l + ... (29) 
4 62a 

The expressions (|28|) and (f29|) do not contain a term 
proportional to the area of the shell, which is a sign of 
zero coefficient of surface tension in the Casimir prob- 
lem. The only geometric contribution present is the Kac 
term with K.( £M ^ = 1/4. In contrast to its scalar coun- 
terpart IC^l = -1/24 (see Eqs.(QI]) and (US)) the elec- 
tromagnetic Kac number is positive which means that 
a conductive boundary increases the number of states. 
Our expression for the Kac number is a special case of 
the remarkable prediction due to Balian and Duplantier 
15] who established that in the high-temperature limit 
the energy of a conductive shell is proportional to the 
temperature with a coefficient BT> given by 

BV = — l — [ daiinl + in 2 2 + 2n x n 2 ) - n (30) 

1287T J 

This was interpreted as the number of additional modes 
created by the shell. It depends on the topology of the 
surface through its genus n, and on its local curvatures 
«i,2- For the sphere this gives BT> = 1/4, in agreement 



with our results (|28[) and ([29]) . It implies a T/4 contri- 
bution to the energy at finite temperature T. 

The Casimir energy is the mode sum for the case 
F(q) = q/2. Its value 3/64a is determined by the deriva- 
tive term of (|28[) . The 3 /64a universal answer agrees 
with the evaluation given by Nesterenko and Pirozhenko 
[22I ] , who used the zeta function regularization procedure. 
Even though this represents the bulk of the Casimir en- 
ergy for a conductive spherical shell [Hf, there are 1/a 
corrections to this result due to higher order terms in the 
asymptotic expression (j2~4"|) . These corrections do not re- 
quire a cutoff and have been calculated (32| . 

IV. SUMMARY 

We have shown how to calculate the generic mode sum 
for three cases: a scalar field in one dimension perturbed 
by a Dirichlet boundary; a scalar field in three dimen- 
sions, perturbed by the introduction of a spherical shell; 
and the electromagnetic field in three dimensions, per- 
turbed by the introduction of a conducting shell. We 
have shown that these will in general contain contribu- 
tions of geometric origin (the Weyl terms) which require 
a cutoff and have physical meaning. We have shown how 
to extract the Kac number and the Casimir term when 
it exists. 
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